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DETERMINATION OF STATISTICS FOR ANY ROTATION OF AXES 
OF A BIVARIATE NORMAL ELLIPTICAL DISTRIBUTION 

SUMMARY 


A method is presented for the determination of the estimates of 
means 5 variances, covariances, and correlations for any rotation of axes 
of a bivariate normal elliptical distribution. This problem arises in 
connection with vehicle launches at Cape Kennedy, Florida, when winds 
normal and parallel to the vehicle flight path are launch constraints. 

A vector wind data sample is resolved into wind components in a recfcangu- 
lar coordinate system, and the bivariate normal elliptical distribution 
theory is applied to the data sample of component winds. However, 
theoretical inferences are required for any arbitrary firing azimuth of 
the launch vehicle. This requires a rotation of axes of the bivariate 
normal surface and the computation of new statistics after rotation. 


I . INTRODUCTION 


In 1846, Bravais [1] presented the theory of normal frequency dis- 
tributions in two and more variables. Most modern s tatis tical textbooks 
include some discussion of such distributions. Outstanding contributions 
are notable in the works of Bertrand [2], Czuber [3], Pearson [4], 

Kluyver [5], and Student [6] „ 

Further developments and some practical applications of these dis- 
tributions with particular regard to geophysical problems are presented 
by Crutcher and Baer [7], Crutcher and Moses [8], and Groenewoud, Hoag 1 in, 
Vitalis , and Crutcher [9]. Anderson [10] and Hald [11] treat the subject 
in general form. This report presents a method for the determination of 
the estimates of means, variances, covariances, and correlations for any 
rotation of axes of a bivariate normal elliptical distribution. Although 
the procedure can be extended to the case for more than two variables, 
this case is not discussed here. 

This problem arises in the determination of probabilities of success- 
ful launches at Cape Kennedy, Florida. For example, winds exceeding some 
given value normal and parallel to the flight path may be defined as 
launch constraints. Those winds normal to the flight path are designated 
as right and left cross winds , while those parallel to the flight path 



are called head and tail winds, A vector wind data sample is resolved 
into wind components in a rectangular coordinate system and the bivariate 
normal elliptical distribution theory is applied to the data sample of 
component winds. However, theoretical inferences are required for any 
arbitrary firing azimuth of the launch vehicle. Making such inferences 
requires a rotation of the axes of the bivariate normal surface and the 
computation of new statistics after rotation. 


II. THE BIVARIATE NORMAL ELLIPTICAL DISTRIBUTION 


Let x x and x 2 be normally distributed variables with parameters 
(i± 3 a ±) and (| 2 ,c r 2 ). £1 and | 2 are the respective means , 

while Q x and cr 2 are the respective standard deviations. Let p be the 
correlation coefficient, which is a measure of the dependence between 
x x and x 2 . Now, the equation of the two-dimensional normal density 
function takes the form 


p( x i* x 2 ) « (2 it tf x a 2 [l~p 2 ] 1 / 2 )~ 1 exp i~(2[l-p 2 ]) 



2p 


&r) - (=£*) ] }. 


Given the variables x x and x 2 , we calculate the estimates of the 
bivariate normal parameters. Denote the estimates as X x , X 2 , S x , 
and r XlX2 , A special case of the problem described here is to 
determine the estimates of the parameters where a rotation of the axes 
through an angle QL reduces the correlation of components to zero. This 
case is discussed in Crutcher and Baer [7], and Groenewoud, Hoaglin, 
Vitalis, and Crutcher [9]. From the statistics in the (x x ,x 2 ) space, 
we must determine the statistics for any rotation of the axes through 
an arbitrary angle a in the following way. 

Let (y x ,y 2 ) be the new space after rotation with orthogonal axes 

yi and y 2 . 


Y x = X x cos QL 4* X 2 sin a (i) 

Y 2 = X 2 cos a - X x sin a, (2) 
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where a is the angle of rotation, and x x and x 2 are the original axes 
along which the values of X x and X 2 are measured , Denote the sample 
variances of (x l5 x 2 ) as S 2 and S 2 , respectively. The covariance of 
(Xi,X 2 ) is 


S = r S S 

x l x 2 x l x 2 X 1 x 2 


For the special case when S 

x l x 2 


covariance matrix is 


or r 


x l x 2 


is equal to zero, the variance- 





— 

s 2 s 


1 

0 

Xi XiX 2 

= A 



S S 2 


0 

1 

X 2 Xi x 2 






S 

X 2 X 1 


s 

x l x 2 



Setting the determinant of this matrix equal to zero and solving 
for the two roots of the quadratic equation is essentially the rotation 
of axes through an angle OC to where components along the new axes are 
uncorrelated. Thus, the roots Ax and A 2 are the variances along the 
major and minor axes, respectively, of the bivariate distribution and 
the covariances (correlations) are zero. From these data, variances 
and covariances (correlations) along any set of new axes can be deter™ 
mined by rotating the axes from this basic position. Then the matrix 
will be 
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2 

s 

y 1 

YlY 2 


s 2 

Y 2 Yl 

y 2 


where the new space is (71,72) and y x and 72 are the axes. 

The statistics in the (71,72) space ma7 be derived as follows: Let 

Yi, Yp, , S 2 , S , r be the means, variances, covariance, and 
" 5 * Yi 7s 7 i72 7i72 

correlation coefficient in the (71,72) space. Taking mathematical 
expectations of equations ( 1 ) and ( 2 ) , we have for the new means 


Y x = X x cos QC + Xp 
Yp = Xp cos OC ” X x 


sin a ( 3 ) 

sin a, ( 4 ) 


Subtracting equation ( 3 ) from equation ( 1 ), we have 


(¥1 - Y x ) = (X x - X x ) cos a + (Xp - Xp) sin a. 


( 5 ) 


and subtracting equation ( 4 ) from equation ( 2 ), we have 

(Y p “ Yp) — (Xp — Xp) cos OC “ (X x X x ) s in Gb ■ ( 6 ) 

Let 

x i “ (Xi - X x ) 

xp = (Xp - Xp) 

71 - Ofi - Y x ) 

72 = (Y 2 - Yp) 


be the deviations from the respective means. Now, equations ( 5 ) and 

( 6 ) become 
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Yi ~ x i cos a -f x 2 sin O' 


(7) 


y 2 = cos a - x x sin a. (8) 


Squaring equations (7) and (8) gives 

2 2 O 2 «~a 

yi = x x cos a + x 2 sin a + 2x x x 2 cos a sin a (9) 

yl = xf cos 2 a + xf sin 2 a - 2xiX 2 cos a sin a. (10) 

Now, taking expected values of equations (9) and (10), we have 
for the new variances in the (y x ,y 2 ) space 

S 2 = S 2 cos 2 Q! + S 2 sin 2 a + 2S cos a sin a (11) 

yi x X x 2 X X X 2 

and 


S 2 = S 2 cos 2 a + S 2 sin 2 a - 2S cos a sin a. (12) 

y 2 x 2 x x x x x 2 ' 


To derive the covariance, Sy^y^, in the (y x ,y 2 ) space, we multiply 
equation (7) by equation (8) and obtain 

yxy 2 = x x x 2 (cos 2 a - sin 2 a) + cos a sin a(x| - x 2 ). (13) 

Taking the mathematical expectation of equation (13) gives the 
covariance in the (y x , y 2 ) space 

S = S (cos 2 a- sin 2 a) + cos a sin a( S 2 - S 2 ) , (14) 

y x y 2 x x x 2 v x 2 


Finally, the correlation coefficient in the (y x ,y 2 ) space is defined 
as 


r 

yxys 


s yxy 2 

s s 
yx y2 


(15) 
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Equations (3), (4), (11), (12), (14), and (15) give the required 
estimates of the statistics in the (y l5 y 2 ) space. 

A computer program whose input is the statistics in the (x l5 x 2 ) 
space and any arbitrary angle of rotation of the bivariate normal sur- 
face is available upon request. The output of the program includes 
the required statistics in the (yi,ys) space. 
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